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1 Introduction 


Streamwise vortices are known to be an important ingredient of the transition process in 
boundary layer flows over both curved and flat walls. In the case of flat walls the streamwise 
vortices are an initially passive product of wave interactions, see for example Hall and Smith 
(1991), whilst for curved walls the vortex can be the primary instability of the undisturbed flow. 
In recent years there has been much interest in the breakdown process of such vortices. In the 
context of Gortler vortices the experiments of Swearingen and Blackwelder (1987) showed that 
the initially steady vortex structure caused by wall curvature undergoes a secondary bifurcation 
to an unsteady three-dimensional flow. This secondary instability can be of either the sinuous or 
varicose type. In flat plate boundary layers the secondary instability is invariably found to be of 
the varicose type. 

The first attempt to understand the breakdown process for Gortler vortices was given by Hall 
and Seddougui (1989). In that calculation the small wavelength Gortler vortex flow investigated 
by Hall and Lakin (1988) was localized in the shear layers trapping the region of vortex activity. 
Analytical progress for this case is made possible by the assumption of small vortex wavelength 
and Hall and Seddougui (1989) showed that wavy (i.e. sinuous) modes of instability are possible 
and are localized near the top and bottom of the vortices. No mode of the varicose type was 
found to be unstable. If such a mode exists at small vortex wavelengths then it is presumably not 
trapped in the region where the vortex activity decays to zero. 

In order to understand the origin of this secondary instability mechanism, Hall and Horseman 
(1991) investigated the inviscid linear instability of a Gortler vortex. The instability analysis for 
such a flow is greatly simplified by the fact that the streamwise velocity component of the flow in 
the presence of a Gortler vortex remains an order of magnitude larger than the normal and spanwise 
velocity components. This means that an inviscid disturbance to the vortex flow satisfies a two- 
dimensional form of the Rayleigh equation dependent only on the streamwise velocity component 
of the vortex. Note here that the lengthscale and timescale of the inviscid mode are relatively 
short so that nonparallel effects are formally negligible in the leading order inviscid instability 
analysis. The modified two-dimensional Rayleigh equation found by Hall and Horseman (1991) 
was simultaneously found in the context of vortex wave interaction theory by Hall and Smith 
(1991). The numerical solution of the eigenvalue problem associated with the modified Rayleigh 
equation was discussed by Hall and Horseman (1991) who found that both sinuous and varicose 
modes become unstable as the vortex develops in the streamwise direction. The question of which 
mode is the most unstable is a function of the history of the vortex, its wavelength and the Gortler 
number. However, Hall and Horseman (1991) were able to obtain quantitative agreement with the 
experimental measurements of Swearingen and Blackwelder (1987). Subsequently the instability 
problem was examined again by Li and Malik (1994) who found that Hall and Horseman (1991) 
had missed the most unstable varicose mode. Note that this latter mode is in fact the second most 
unstable overall and that otherwise the results of Li and Malik (1994) were consistent with those 
of Hall and Horseman (1991). 

Related investigations of the problem discussed above have been given by Liu and Domardzki 
(1993) and Yu and Liu (1994). In the former paper a Direct Navier-Stokes (DNS) simulation of 
transition in the Gortler problem was carried out for parallel flows and the results obtained are 
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in broad agreement with those of Hall and Horseman (1991). Yu and Liu (1994) reconsidered 
the instability problem of Hall and Horseman (1991) but retained some viscous terms in their 
approximation. Yu and Liu (1994) are critical of Hall and Horseman (1991) because of their 
neglect of viscosity. However, this criticism is perhaps surprising since viscous effects are certainly 
a second order effect and in any case the equations solved by Yu and Liu (1994) retain only some of 
the second order effects. In effect the criticism of Yu and Liu (1994) is equivalent arguing that the 
Orr-Sommerfeld equation rather than the Rayleigh equation should be used to describe inviscid 
instabilities of parallel flows. 

In this paper we shall describe the evolution of the inviscid mode found by Hall and Horseman 
(1991). This will be done using viscous-critical-layer and diffusion-layer theories in the context 
of a weakly nonlinear instability theory. In particular, we shall consider the evolution of a mode 
near the critical streamwise location where the vortex structure has developed sufficiently to (first) 
render the now three-dimensional boundary-layer flow unstable to inviscid modes (note that, for 
example, incompressible two-dimensional Blasius-boundary-layer flow does not support inviscid 
instability). At such a location, the flow is marginally unstable and we can consider the evolution 
of the most dangerous (important) mode. Related work has been carried out by Wu (1993) and 
Smith, Brown & Brown (1993). 

Although the analysis we give is for vortex flows generated by wall curvature, it is valid for 
any flow where one of the velocity components depends on two spatial variables and is larger than 
the other two components. For such more general flows the periodicity in the spanwise direction 
which we assume in this paper must be replaced by an appropriate condition in order to derive 
the required solvability condition. The procedure adopted in the rest of this paper is as follows: 
In Section 2 we formulate the problem to be considered; in Section 3 we determine the outer 
solution of the perturbation equations and the required form of the solvability condition, while 
the critical layer is discussed in Section 4. At the edge of the critical layer a diffusion layer is 
required in order to account for the mean flow correction; this layer is discussed in Section 5 and 
the evolution equation is derived; the solution of the evolution equation and some conclusions are 
given in Section 6. 


2 Formulation 

The nondimensionalised Navier-Stokes equations for an incompressible flow may be written 
in the form 


U X + Vy + W Z} 

u t + uu x + vu y + wu z = -p x + + u yy + u zz ), 

V t + UV X + vv y + wv z = -p y + Re^iyxx + v yy + v zz ), 
w t + uw x + vw y + ww z = -p z + Re~ 1 (w xx + w yy + w zz ) 
where Re is the Reynolds number. 


(2.1a — d) 
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We define new variables Y, Z by writing 

y = R e ~ 1/2 Y . , 2 = Re~ 1/2 Z } (2.2 a, 6) 

and assume a three-dimensional boundary-layer flow (i.e. Blasius flow plus streamwise vortex) of 
the form: 

(u,v,w,p) = (u(x,Y, Z), i?e _1/,2 n(x, Y, Z), Re^^w^x, Y, Z), Poo + Re~ 1 p(x , Y, Z)). 

This flow can be generated by several mechanisms but the most obvious one is wall curvature. In 
that case x,y in (2.1) measure distance along and normal to the wall and curvature terms must 
be inserted into those equations. We then find that the flow is determined by 

u x + v Y + W Z , 

UU X + VUy + WUz = UYY + U ZZ , 

i„ 2 

Vt + UV X + VVy + wv z + -Gu = -py + Vyy + V ZZ , 

U) t + UW X + VWy + WW Z = —p Z + Wyy + W ZZ ] (2.3a - e) 

which must be solved in conjunction with suitable boundary conditions such as ‘no-slip’ at the 
wall (solid boundary) and with the solution tending towards a uniform flow at infinity (far from 
the solid boundary). The parameter G appearing in (2.3) is the Gortler number. The extra 
curvature terms to be inserted in (2.1) play no role in the nonlinear inviscid instability problem 
to be investigated here so for convenience we choose not to insert them. Note also that, from Hall 
and Bennett (1986), the Gortler equations are perhaps more conveniently derived by starting with 
(2.1) and then making an appropriate Prandtl transformation. Again if this route is followed the 
extra terms to be inserted into (2.1) play no role in the following analysis. Thus in the present 
investigation it is sufficient for us to use (2.1) to describe the nonlinear state perturbed around the 
incoming vortex flow given above. Whether the vortex held is generated by curvature, turbulence 
screens upstream, or a localized bump at the wall is irrelevant in the following discussion. 

In general the solution of (2.3) is a numerical task, see Hall (1988), and it turns out that non- 
linear effects stabilize the growth of Gortler vortices; however when a large spanwise wavenumber 
assumption is made, Hall & Lakin (1988) demonstrated that much analytical progress can be 
made towards the solution of these equations. 

Let us first recap the linear inviscid stability problem for this three-dimensional boundary- 
layer flow. In the neighbourhood of a point x 0 the flow is perturbed by a small inviscid disturbance 
proportional to 

E = exp [i(aX — OT)], 

where 

x = x 0 + Re~ 1 ^ 2 X and t = i?e _1/,2 T, (2.4a — c) 

and a, fl are the (streamwise) wavenumber and frequency of the linear, inviscid secondary insta- 
bility. The expansions for the velocities and pressure are 

(a, v, w,p) = (a, i?e -1/,2 n, i?e _1/,2 rc, p^ + i?e _1 p(x, Y, Z)) + e ((a, v, w,p)E + c.c.) , (2-5) 
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where ‘c.c.’ represents complex conjugate; barred quantities correspond to the three-dimensional 
boundary-layer flow; and the disturbance quantities u } v } w } p are, in particular, functions of x,Y 
and Z (but not X or T). 

After a little algebra we find that the pressure perturbation p satisfies the modified Rayleigh 
pressure equation 

d_ 
dY 

with boundary conditions 

p Y (y = 0) = 0, p(Y — >■ oo) = 0, (2.66, c) 

where c = ft /a. This equation was derived by Hall and Smith (1991) who were concerned with 
vortex-wave interactions and by Hall and Horseman (1991) in the context of secondary instabilities 
of Gortler vortices. The eigenvalue problem for c = c(u, x, a) (temporal stability problem) associ- 
ated with the partial differential system for p was first solved by Hall and Horseman (1991). Here 
we shall consider the more appropriate spatial instability problem in the presence of nonlinear 
effects. 

Note that c is not a function of Z; if we consider neutral disturbances (those having c entirely 
real), equation (2.6a) is singular at Y = Yc = f(x,Z ), (say) where u = c. Thus for three- 
dimensional boundary-layers, the critical layer is ‘wavy’ in the sense that the location of the 
critical level (where the equation is singular) is a function of spanwise location Z . Note further 
that the neutral value of c, and hence Yc = f(x,Z ), are not known in advance of a numerical 
solution of the eigenvalue-problem (2.6a-c). 

In Figure 1 we show some results for the most dangerous odd mode of instability for a Gortler 
vortex flow. The shape of the curve is similar to that which would be obtained for an inflectional 
uni-directional velocity held by solving the standard Rayleigh equation. In Figure 2 we show the 
normal and streamwise perturbation velocity components for an unstable mode. Note that the 
mode is concentrated towards the edge of the boundary layer. Since the mode is not neutral it does 
not exhibit a critical layer behaviour. In fact it is a difficult numerical procedure to solve (2.6) in 
the neutral case and results are not yet available. Now let us consider the nonlinear problem and 
as a starting point we discuss the outer region away from the critical layer. 


Py 


[u — c 


d_ 
+ dZ 


Pz 


[u — c 


a 2 p 
(u — c) 


= 0 , 


(2.6a) 


3 The outer solution for weakly nonlinear inviscid modes 


In order to derive the desired nonlinear evolution equation for the amplitude of an inviscid 
disturbance mode, it is necessary to split the three-dimensional boundary-layer how into distinct 
regions (layers), each corresponding to different dominant physical effects locally governing the 
inviscid disturbance; see Figure 3. Let us hrst consider the how solutions in regions Ia,b, away from 
the critical layer, diffusion layer and solid boundaries (the wall) i.e. the bulk of the boundary- 
layer. Our aim is to derive the solvability condition associated with an inhomogeneous form of 
(2.6a-c). This solvability condition, together with expressions for certain ‘jumps’ in how velocities 
in crossing the critical layer, (to be derived in the later subsections) will yield an evolution equation 
for the spatial evolution of the disturbance amplitude. 


4 



Before proceeding any further, let us consider the various streamwise lengthscales which are 
present in the analysis. In fact, we employ multiple scales in the streamwise direction x: 

x = x 0 + fJ-x, x = x 0 + Re~ x ^ 2 /r~ x X } x = x 0 + Re~ x / 2 X. (3-1) 

Here x 0 now denotes the neutral ^-station (in fact, x 0 is a point of marginal instability as the 
flow first becomes unstable to linear inviscid disturbances here); fix is a small change in £ -location 
from x 0 (assuming fix > 0, then the flow is now slightly unstable: hence the weakly nonlinear 
analysis to follow); Re~ x A p~ x X is the wave-amplitude modulation scale; and Re~ x / 2 X is the scale 
on which the waves oscillate. Note that we are considering the purely spatial stability problem 
and thus there are no slow time effects. 

We consider ‘fixed-frequency’ disturbances so that 0 is fixed; thus the non-neutrality is entirely 
due to the the change in £ -location gx. Note that the developing three-dimensional boundary- 
layer, and hence the neutral streamwise-wavenumber a, are functions of x. In addition it should be 
noted that, if we were considering a three-dimensional shear-layer rather than a three-dimensional 
boundary-layer, we would not consider 0 fixed. The non-neutrality would then be entirely due 
to the 0 perturbation (i.e. a Strouhal number perturbation). 

Since the normal and spanwise velocity components of the unperturbed flow are relatively 
small only the streamwise component u of the three-dimensional boundary-layer flow enters the 
problem to orders of concern. We write 

u = u 0 -\- /rui + ..., u 0 (Y } Z) = u(xo, Y, Z), ui(Y, Z) = u 0x (x 0 , Y, Z)x } • • • . (3.2a) 

Further, we only need to consider two terms of the fundamental (proportional to E ± ) and write 

( u,v,w,p ) = (ui,vi,wi,pi) H b p{u 2 ,V 2 ,w 2 ,p 2 ) H • (3.26) 


(i) The leading fundamental term 

The leading order problem for (hi, hi, hq,pi) is a partial differential system in Y, Z and its 
solution may be written in the form 

(hi, hi, wi,pi) = A(X)(v^\ v[ x \ w[ x \p^)E + c.c., (3.3a) 

where A is an amplitude function and h^, hj 1 ^, rhj 1 ^ and satisfy 

■ W) , W) , W) n ■ (- w(i) , - NT , - W) ■ dB 
iau\ + v]y + w\z = U, ia[u 0 — c)u\ + Uqyv \ + Uqzw\ = — iap\ 

ia(u 0 — c)^! 1 ^ = — Piy, ia(u 0 — c)^^ = —p{z- (3.3 b — e) 

These lead to the eigenvalue problem (2.6) (with replacing p) to determine the neutral solution 
at x = xq. We fix p^ by the normalization: 

= poo(Z)e~ aY as Y — >■ oo, (3.3/) 
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where p^Z) is prescribed at some value of Z. 

Following Hall and Smith (1991) let us now consider the behaviours of ..., near the 
critical layer at Y = Y c = f(Z). Defining s = Y — Y c (so that dz — > dz ~ fzd s ) } we write 

uq — c A(x, Z)s Y X 2 (x } Z)s 2 /2 T Ag(x, Z^js 2 /6 -)-•••, 


U\ — U\ q(x ^ Z) Y ^ 11(^7 Z)s Y 2 (x, Z)s / 2 Y * * * 
The method of Frobenius gives in the form 

ph = 4>i + bMZ)h, 


( 3 . 4 ) 


(3.5a) 


where 


<?h — <f>io(Z) + ^ , ii(-^)- 5 + <^ 12 (Z)-s 2 Y (f>i 3 L(Z).s 3 In 5 + O.s 3 + • • • , (f > 2 — s 3 Y • • • • (3.56) 

Note that <p w (Z) and bi± have to be determined numerically (the ± corresponding to above 
and below critical level respectively); while, in particular, 

A</>n = fz<f>ioz, 

2AA 2 (f>i 2 = XA.(f) WZ z ~ (2A Z + Xf z fzz)4 , ioz — a 2 XA<f> 10 

and 

6AA^i 3i = XXf Z (f>i 2 Z + 2(3A 2 A — ifzXz Y 2 Xf Z z )</>i 2 — 2A (f>nzz + 2(2X Z — X 2 fz)<t > nz 

Y{2a 2 X — 2f z X 2Z + X 2 f Z z)<t>n ~ X 2 <f> w zz + 2A 2 ^</>i 0 z + at 2 X 2 cf> w , (3.6a — c) 

where we have defined 

A = 1 Yfl (3.7) 

Here, following Hall and Smith (1991), we take <f>\ 3 L = 0. If it were non-zero, at 0(Re~ 1 ^ 6 ) 
a jump in v across the sole critical layer would be induced and as we are not looking at upper- 
branch scalings there would be no wall layer effects to counteract this jump; hence this jump and 
the coefficient of the logarithm-term must be zero. 

The velocities and have the following properties near the critical layer: 

= d 1 zs ~ 1 Y • • • , Vpl = — fzW\^ = d 2 Y • • • , = — iadis~ x + • • • , 

where 

di = —<f>wz/ a 2 AA and d 2 = (— 2A^i 2 + fz<t>iiz)/i(* A. (3.8a — e) 
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(ii) The largest forced fundamental 


Let us now consider the largest forced term of the fundamental (due to non-neutrality effects), 
denoted by (u 9 \ w 9 \ p' 9 )- It is found that 


9 PzY 9 P { 2Z a2 P2 } _ R 2 

dY ( u-c ) 2 + dZ ( u-c ) 2 (u — c) 2 (u — c) 2 ’ 


(3.9a) 


with boundary conditions 

P 2 y( y = 0) = 0, p { 2 ] {Y -> 00 ) = 0, (3.96, c) 

and 

R2 = R ( 2 ) ^ + R ( 2 ) A, (3.9 d) 

oX 

where 

R ( 2 ] = —r. — -APiyy + Pizz ~ « 2 M 1} ) - 2 iap? , 

a[u 0 — c) 

R? = ~ ( - Ul A Piyy + Pizz ~ a2 Pi ] ) ~ J— — r(“i yPiy + u\ zPiz) ■ ( 3 - 9e , /) 

_(u 0 -c) ( UQ-C ) 

Note that the equation for p' 9 is an inhomogeneous form of that for as mentioned previ- 
ously, we must determine a solvability condition for this equation, to ensure that it has a solution. 
This solvability condition will (indirectly) lead to the desired evolution equation. To derive the 
solvability condition we will essentially follow the conventional method of multiplying the inhomo- 
geneous equation by the adjoint of the homogeneous equation and then integrate over the range 
0 < Y < 00 , 0 < Z < 27T (the latter corresponding to a complete period). Note that we are 
assuming that the disturbance has the same period as the basic flow. The latter assumption is 
justified by the fact that Hall and Horseman (1991) were unable to find any subharmonic distur- 
bances of the modified Rayleigh pressure equation. As in all critical layer problems, special care 
must be taken to deal with the singular nature of the equations (i.e. the critical layers); however 
extra care is necessary for the current problem due to the ‘waviness’ (Z-dependence of Y c = f(Z )) 
of the critical layer. The fact that additional care is required can be explained as follows: note 
that in (Y, Z)-coordinates (i) for any fixed Z-value there is a unique critical-Y value, but (ii), for 
a fixed Y-value there is either no critical Z-values or many. 

Thus we look for co-ordinates which describe the critical-level more suitably; note that the 
Prandtl-transformation (5 = Y — f(Z )) does ‘level’ out the critical-layer but, at the same time, 
leaves the solid boundary ‘wavy’ (i.e. it is described by Y = —f(Z)). Instead we introduce the 
normal variable 

£ = Y/f(Z); (3.10) 

this transformation ‘flattens’ out the critical-level (it corresponds to £ = 1) and also leaves the 
wall flat (£ = 0). Note that this transformation involves f(Z) and thus could not have been used 
from the outset i.e. it involves information from the solution of the homogeneous problem (1.6a-c). 
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Note that 


d/d Y - (l /f(z))d/d c - d/d z - d/d z + g(z)Cd/dy 

where g = — / //. In these co-ordinates, the homogeneous equation takes the form: 


(3.11) 


( l // 2 


d 


h (1) 

Pn 


it — c 


+ [d/dz + g(d/d c ) 


(d/dz+nHd/90^ 


It — C 


2 A (l) 

crffi ; 

(u — cV 


= 0 , 


= 0) = 0, = oo) = 0; 

whilst the adjoint, q say, satisfies 


(3.12a - c) 


( l // 2 
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it — c 


+ {d/dz - gd/dti) 


(d/dz - gd/d c ()q 


it — c 


a 2 q 


it — c 


= 0, 


<u(£ = 0) = 0, q(£ = oo) = 0. (3.13a - c) 

Therefore in these transformed co-ordinates the pressure equation is not self-adjoint. 

Thus to derive the solvability condition for (3.9a-f) we first transform to (£, Z) variables; 
multiply both sides by the adjoint q (defined above) and then integrate both sides over the range 
0 < C < °°, 0 < Z < 27T, excluding the critical layer. This gives 


p2-7T 


dZ 


1 


^ c y \Cg(qp ( 2z ~ p ( 2 ] qz ) + jziqPx ~ P ( 2 )( id + g 2 d 2 qp ( n ~ £ p 2 1} (£?)«) 


i t=i + 
t=i- 


f 2v qR 2 
Jo Jo I 


[u — c 


7 d^dZ, 


( 3 . 14 ) 


where the bar through the integral represents the finite part (i.e. excluding critical layer effects). 
Near the critical level, q is of the form 

q = qw + <Zn(£ _ 1) + <Zi 2 (£ _ l) 2 + 0(£ — l) 3 ln(£ — 1) + <?i3(£ ~ l) 3 + ' ' ' ? 

whilst i? 2 °\ have the forms 

R { 2 ] = - I)” 1 + r? } + r[ k) (C ~ 1) + r^(£ - l) 2 + • • • , k = 0, 1, 


(the expressions for the coefficients are simply obtained by expanding (3.9e,f) near the critical 
layer, however they are somewhat long and for the sake of brevity are not given here) while 
has the form 

P2 ] = a 2<t>l + b 2 ±<f>2 + P2PI- 

Here the first two terms correspond to solutions of the homogeneous equation, a 2 ,& 2 ± are as yet 
undetermined constants, and subscript PI denotes particular integral. We define a 2 ,& 2 ± so that 
p^pi has no terms in (£ — 1)°, (£ — l) 3 ; i.e. p^pi is properly defined. Then, near the critical layer, 
p^pi has the form 

P 2 .p 1 = fJi(t ~ 1) + fad ~ l) 2 + (JohH ~ l) 3 ln(C — 1) + 0(C — l) 3 + ' ' ' ? 
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where 


Al = /3h’|i + &A, 

oX 


and 


Pzl = T 7 (?i 2 r 3 + qiir ( 0 k) + ?iorf°) - ^( 911 ^ A 2 + q w r {k) A 2 ) + gl °^ 1 (^f - |^), k = 0, 1 


^2 —1 1 111* U 1 1-LU- 1 / ^3 

Thus the solvability condition becomes 

3A q 




qR 2 


2A 2 2A 2 


: d£dZ. 


3A 


( 3 . 15 ) 


/o f 2 A 2 Jo Jo (u — cj 

Near the critical layer, the Prandtl-transformed normal velocity has the form 
V { 2 P = - fzW ( 2 ] = 7T_i((f - 1) _1 + 7T 0 + TT 1L (£ ~ l)ln(£ - 1) + 7Ti±(^ - 1) + • • • (£ > 1) 

where 

3A 3A 

Tr 1L = -——f3 3L and k 1+ - 7Ti_ = ~(b 2+ - b 2 -). (3.16 a - c) 

IQA IQ. A 

Note that for £ < 1, the logarithm ln(£ — 1) in (3.16a) must be replaced by In |£ — 1| — iir; the 
negative sign in front of iir follows from an inspection of the Stokes lines of the operator that 
occurs in the governing equations for the linear problem in the critical-layer. 

Thus, the linear part of the jump that suffers across the critical layer is 


N 1 ) 1 1+ , ■ 

= %1 + ~ + Z7r7r l L] 


( 3 . 17 ) 


an expression for the nonlinear contribution to the jump is calculated in the following two subsec- 
tions before the two expressions are combined to give the evolution equation. 


4 The Critical Layer 


Let us now consider the critical-layer flow (corresponding to region III in Figure 3) with our 
aim to calculate expressions for the jumps in the solvability condition obtained in the previous 
subsection; will shall find that another region, the so-called diffusion layer, also needs to be 
considered and this is the subject of the next subsection. The analysis has similarities to that 
given in Appendix B of Hall & Smith (1991); Appendix A of Brown, Brown, Smith & Timoshin 
(1993); and Wu (1993). Therefore only the essential details of the analysis will be given here. 

The critical layer is centred on F = Y c = f(x,Z ) and has thickness 0(Re~ 1 ^ 6 ) relative to 
the boundary-layer thickness; it is a viscous critical layer. We introduce the critical layer normal 
variable r] ~ 0 ( 1 ) where 

Re 1/2 y = Y = f(Z) + Re~ 1/6 r/ } (4.1a) 
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so that now 


(4.16) 


>■ Re 1 / 6 f z — 

dz dz Jz d v 

and the transformed normal velocity vp is dehned by 

vp = v — fzw. (4.1c) 

The underlying three-dimensional flow has the form 

u = c + Re- 1/6 \(Z)rj + Re- 1/3 \ 2 (Z)r ] 2 /2 + • • • + 

fiu 10 fiRe ^ jiRe ^ ^ 12/2 * * * • (4.2) 

Since we are looking at a viscous critical layer (the flow is marginally unstable) we must assume 
that, 

fi <C i?e -1 / 6 . 

However, depending on the size of //, the above expansion (4.2) will be disordered. 

In the critical layer the perturbation to the three-dimensional underlying base flow has the 
following form: 

(u,v P ,w,p) = e[(i?e 1/6 Hi,I>Pi,i?e 1/6 IH 1 ,A)^ + c.c.] 

+ • • • + e 2 {Re 1/2 U 0 , Re 1/3 V P 0 , i?e 1/2 IH 0 , Re 1/3 P 0 )E° + • • • . (4.3) 

Here only the terms of direct interest have been shown. The second harmonic does not contribute 
to the nonlinear jump and therefore has not been included above. 


(i) The fundamental 

Writing 

U t = A (1) + Re~ 1/6 U[ 2) + • • • , (4.4) 

with similar expansions for Vpi 7 W\ and Pi, and substituting into the governing (Prandtl-trans- 
formed) equations leads to the following system for the leading order fundamental 

iaU? + Vff n + = 0 , 

iaXriUP + XVff + XzrfW^ = + AE/g, 

Pff = 0 

iaXifW^ = -Pg + fzPff + AWg, (4.5a - d) 

which must be solved subject to the condition that solutions match to their outer counterparts 
as i] — > Too. The y-momentum equation of the next order problem (for u[ 2 \ etc) yields, in 
particular, 

P}? = SzPh /A, (4,5e) 


10 



enabling a solution of (4.5a-d) to be determined. The solutions which match with the outer flow 
can be found in the paper by Brown, Brown, Smith & Timoshin (1993); in our notation they are 

A (1) = A(A')*o(Z), 


and 


where 




1 dP [ ^ 

a 2/3 A 2 Q Z J 0 eX P 


icP^rjt — t 3 / 3 dt, 


i dW{ 1] i 
a dZ 4 a 


A z 

A a 


d_ 

dr/ 


W-r. (1) ), 


v>(i) ( ^z , «z\ r/W^ t 1 { X z fzfzz 


a = Xa/ A. 


(4.6) 

(4.7) 


At this stage we also point out that A, a, a, A are all real. 

It is also possible to consider further terms in the fundamental; in fact, by considering the 
terms due to the non-neutrality it is possible to derive the logarithmic jump effect deduced earlier 
from the asymptotic form of the outer solution as the critical layer is approached. As a number of 
different orders have to be considered, and the answers are algebraically messy, we do not present 
that analysis here. 


(ii) The zeroth harmonic 


We now consider the largest zeroth harmonic, or mean-flow correction, (proportional to E ° ) 
due to the nonlinear interactions of the fundamental. Writing 

U 0 = U™ + Re^U^ + • • • , (4.8) 

with similar expansions for Vpo 7 Wo and P 0 , and substituting into the governing (Prandtl-trans- 
formed) equations leads to the following system for the leading order zeroth harmonic 

VJS, + "n/ = o, 

A &ff v - AVpo } - A zrjW^ = S 1} 

Pff = 0, 

f Z P < S - A Pff - 2 f zz w[ 1] w[ lV = 0, 

AH,;.' 1 + fzPff ~ P { S = S 2 , 

where, in particular, 

S 2 = laU^W^ + + w[ x) *w[^ + C.C.. (4.9a - /) 
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As usual in such studies we choose P^ = 0 (any non-zero choice would just correspond to 


different mean-flow). Thus 


from which it follows that 


awQ, = S 2 + 2 j -^j^wp , wp , \ 


fzfzz 


d 1 ) Ad 1 )* 


W r 


(i) 


[±F W + Or , , ~[±F W + OzV 2 /2, 


and 


where 


^ 1) ^24Aa( A 2[± ^ + C ' ] )^ 4 ’ aS 11 °°’ 


{5a z , 7A Z 


F W (X,Z)= ( \<f> 10 \*) z -t-± + 


V 3 a 2A 


no 


tt( 2/3) 2 / 3 T(1/3) 
a 5 / 3 A 5 


|A(A)p 


(4.10a — d) 


and it is important at this stage to note that the X and Z dependences of F w are separable. 

These asymptotes imply that the zeroth harmonic grows on leaving the critical-layer and that 
it is necessary to consider another flow region (the diffusion layer) where diffusion effects can 
counteract this growth. In general such diffusion-layers need to be considered in all nonlinear 
analyses dealing with viscous critical layers (see, for example, Brown & Stewartson, 1978; Wu, 
1993; cf. the buffer-layer of weakly nonlinear vortex-wave theories). 


5 The Diffusion Layer and the Evolution Equation 


This layer (regions IIa,b in Figure 3) is introduced to take care of the growing mean-flow 
corrections as the critical layer is left; it turns out that the nonlinear term in the desired evolution 
equation for A stems entirely from this region. It is necessary to introduce the new scaled normal 
coordinate 

fj = (Y - f(Z))l(6 2 X 1/2 c~ 1/2 ), S 2 = Re~ 1/4 g~ 1/2 < 1; (5.1a, b) 

and so here the mean flow has the form 

u = c + (A/c) 1 / 2 6 2 \rj + • • • . (5.1c) 

The perturbation of this mean flow has the form 

(a, up, w,p) = e [(h 2 1 Ui, Vpi, S 2 1 !l / 7 i, P\)E + c.c.\ + • • • + 

e 2 (S 2 Re 7 / 6 Uo, 6 2 2 Re 2 ' 3 V P0 , 6 2 Re 2 ' 3 W 0 , 6 2 1 P 0 )E° + • • • . (5.2) 

The fundamental has the expansion 

Ui = u[ 1] + • • • + S 3 ui 2) + • • • , h 3 = e 2 h 3 i?e 7 / 6 , (5.3) 
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with similar series for Vpo, Wo and P 0 . Note that the largest fundamental in this region is merely 
a continuation of the outer solution in the locality of the critical layer: 


m 


a) 


®,^i ( Vi (1) ) = (c 1/2 A- 1/2 ^ 


d 2 , —iac 1 ! 2 A x ! 2 d\f] x , (f> w ) A(X) 


(5.4) 


(i) The leading zeroth harmonic 


Writing 

(Uo, Vpo, Wo) = (Ui>\ V^M y ) + • • • , 

we find that the governing equations for the leading zeroth harmonic are 


+ UfJ = 0. 


' a 

ax 


(P_ 

dfj 2 


of + -V% + = o, 


A) 


A X / 2 A; 


(i) _ 


and 


' d_ 
dX 


(P_ 

dfj 2 


Wn (1) = 0. 


(5.5a — c) 


It follows, from matching with the critical-layer solutions, that the boundary conditions at 
fj = ±0 are 


where F± = 
give 






ap = v (A2f±)z ' 


9 ! i'g 


= -F. 


r\ ~cy -*- zt.Zv} r\ ~ 

a /)" 1 of) 

A(± F w + C)/c. These equations can be solved using a Fourier transform method to 


= 


Uoll = ±( 1 / 2 ctt 1/2 )?7 ^ dtt~ 1/2 exp(-?f /4t) (A [F±(X - t)] z 

+ (l/4)^r 1 A z F ± (A -t)- (3/2)A z F ± (X - tj) , 

VpS = ±7T“ 1/2 r 1/2 [F±(A - t)] z exp(-?f/4t)dt, 

IFg 1 ^ = ^(c/pA) -1 / 2 f t~ 1 / 2 F±(X — t) exp(— fj 2 /4:t)dt. (5.6 a — c) 

Jo 

As the meanflow correction is larger in the diffusion layer than the critical layer, the leading 
order nonlinear jump will result from the interaction of the meanflow and the largest fundamental 
in this (the diffusion) layer. 
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(ii) The largest forced fundamental 


The governing equations for the largest fundamental due to non-neutrality effects are 

t a(f> + (c/AJ'Cyyl + Wg = 0. 

ta(A/c) 1/2 X}Ug + A i'g + (A/c) l/2 \ z ,jWg + taPg = 

-uf'i'fx - ( c/a f 2 vgug! - wgu^. 


P$ = 0, 

.Q(A/ C ) 1 /r\,)iif ' + Lpm = -ty'iu'i 1 . 

These can be combined to give an equation for VjJj'. 


= (c/A)'/ 2 \-'A(X) 


d_ 

df} 


-Zictdrf-'iU® ~ -ft!™) + d 2 U$ 


substituting for Uq 1 ^ from equation (5.6) gives the required nonlinear jump 

A 1 / 2 . , ~ s ( . iadx , „ X 


ty( 2 ) 

v Pifj 


A(V\ I J "^“'1 (Pi \ A 1 [°° \ 2 ( F W (X Z) 

—A(X)U-—(dz- T )A\J o A ^ 


(5.7a — d ) 


(5.7e) 


dC- (5-8) 


(ill) The evolution equation 


Noting that F w is separable in X and Z } we introduce F WW (Z ) such that 


F W (X,Z) = F WW (Z) \A(X)\ 2 . 


(5.9 a) 


Equation (5.8) may then be written in the form 



OO 

— oo 


A 1 / 2 

XFA 


d 2 



N)a) A 2 A( .Y) j“ \A(X - C)| 2 dC- 


(5.96) 


Thus the hnal evolution equation can be obtained from (3.15) and (5.9b). If the constants 
7i , 72 are defined by 


7i 


f 2v I [°° qR { 2 ] _ SiwqwA^ 

Jo ( u-c ) 2 * A 2 / 2 



qRj ] 


[u — c 


-d£ 


3iTrq w X/3 


(!) 
3 L 


dZ , 
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72 


27T 


0 


iaqio 

TfP 


d 2 


zcxd\ 

~ 


(dz 


H a2 


F WW (Z) 

A 






then the required evolution equation is 


3i7rgioA/?.^ 

A 2 / 2 


c£Z, 


(5.10a, 6 ) 


W4 „ /-oo 

T- + 7i(*o)A = a - .i'd i. li A ! / \A(X - C)| 2 dC, 
dX Jo 


A^O 


as 


X -> - 


oo. 


( 5 . 11 ) 


Here the parameter 


e 2 i?e 2 / 3 


( 5 . 12 ) 


measures the ‘competitiveness’ of the nonlinear term, in the evolution equation, relative to the 
term corresponding to linear effects. For <r <C 1, linear effects will dominate the evolution and 
the disturbance will continue to grow after x = x 0 . Here we concern ourselves with the regime 
a ~ 0 ( 1 ), where linear and nonlinear effects have equal influence on the evolution of the distur- 
bance. A more general form of (5.11) can be derived which contains an additional ‘linear’ term, 
proportional to XA, to account for non-parallelism effects which are important for certain ranges 
of the parameters (e, //, Re) cf Hall & Smith (1984); Smith, Brown & Brown (1993). 

Thus we see that the Z dependence of the problem has been removed by the application of 
the solvability condition so that the Z dependence of the nonlinear problem is nonlocal as was 
found to be the case in the linear regime. Further, we see that our weakly nonlinear analysis has 
lead to a cubic nonlinearity; however rather than appearing as a polynomial ( eg. as A|A | 2 if the 
evolution was described by the Stuart- Watson method), the nonlinear term is a convolution. The 
evolution equation is an integro-differential equation which depends on the entire history of the 
disturbance. Such evolution equations were first derived/proposed by Hickernell (1984). In fact 
Wu (1993) and Smith, Brown & Brown (1993) have derived essentially the same equation in their 
studies of boundary-layer transition; however, as we have considered a fully three-dimensional 
boundary-layer, our coefficients 71,72 are far more complicated. Similar equations have also 
appeared in other recent papers eg. Smith & Walton (1989), Smith & Blennerhassett (1992), Wu, 
Lee & Cowley (1993) and Blackaby (1994). 


6 Solution of the evolution equation and conclusions 


The nature of the solution of (5.11) depends crucially on the sign of the real part of j 2 
(defined by equation (5.10b). Since the disturbance under investigation becomes unstable as it 
moves downstream we know that the real part of 71 is positive. Without calculating the solution 
of the neutral leading order eigenfunction and adjoint problems we cannot say what is the sign of 
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the real part of y 2 . We shall therefore discuss both possibilities and use results from experiments 
to suggest the most likely scenario. Firstly we note that a suitably rescaled version of (5.10) takes 
the form: 


dH - r°° 

— t = B±B{X) B(X-C) dC, 0 as X -> - 00 . (6.1) 

dX Jo 

Here B is real, positive and proportional to \A\ 2 and the ± signs correspond respectively to the 
cases when the real part of y 2 is positive and negative respectively. 

The solution of (6.1) which has the required upstream behaviour is given by 


B 


2e* 

(1 =p e x ) 2 ' 


( 6 . 2 ) 


If we take the negative sign in (6.2), corresponding to the positive sign in equation (6.1), we 
see that a singularity develops after a finite distance. If the positive sign is taken in (6.2), B grows 
as X increases from —00 until it reaches a maximum and then decays exponentially to zero. At 
first sight this seems a rather curious fate for a disturbance which was initially unstable on the 
basis of linear theory. However the integral term in (6.1) can be interpreted as the effect on the 
growth rate by the mean flow corrected by the upstream development of the instability. Thus the 
mean flow modification which occurs in the early stages of the growth of the disturbance adjusts 
the meanflow so as to make it linearly stable further downstream. 

Experimental observations certainly suggest that the mode identified by Hall and Horseman 
(1991) continues to grow after it first becomes unstable. Here we are assuming that the mode of 
instability discussed by Hall and Horseman (1991) is responsible for the experimentally observed 
onset of three-dimensionality in the Gortler problem. The closeness of the theoretically predicted 
most unstable wavenumber and frequency with those measured by Swearingen and Blackwelder 
(1987) gives some backing for that assumption. The experiments therefore suggest that the nega- 
tive sign is appropriate in (6.2). However without numerical solutions of the neutral eigenfunction 
and its adjoint we cannot confirm that assertion. 

We have carried out a viscous critical layer analysis for a marginally unstable inviscid distur- 
bance to a flow containing a streamwise vortex structure. The vortex structure could be the result 
of a centrifugal instability, wave interactions or other mechanisms. In fact our analysis is valid 
for any flow where one of the velocity components depends on two spatial variables and is larger 
than the the other two components. For such more general flows, the periodicity in the spanwise 
direction, which we assume in this paper, must be replaced by an appropriate condition in order 
to derive the required solvability condition. 

Our analysis is similar to that of Wu (1993) and Smith, Brown & Brown (1993); however 
our analysis is complicated due to the fact that the nonlinear vortex state has rendered the 
boundary-layer flow three-dimensional. As a consequence, our disturbances have a general Z- 
dependence whereas Wu (1993) was able to consider separate harmonics in Z and derives coupled 
amplitude equations. In the problem considered by Smith, Brown & Brown (1993), the initial 
boundary-layer flow is two-dimensional; all subsequent vortex activity (three-dimensionality of 
the boundary-layer) is due to the relatively large mean-flow corrections induced in the diffusion 
(buffer) layers via nonlinear-interaction effects. 
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In this paper we have described the evolution of the inviscid modes found by Hall and Horseman 
(1991); this has done using viscous-critical-layer and diffusion-layer theories in the context of a 
weakly nonlinear instability theory. In particular, we have considered the evolution of a mode 
near the critical streamwise location where the vortex structure has developed sufficiently to 
(first) render the now three-dimensional boundary-layer flow unstable to inviscid modes. At 
such a location, the flow is marginally unstable and we can consider the evolution of the most 
dangerous (important) mode. We note that our theory is not directly applicable to modes excited 
at streamwise locations where the flow supports a band of unstable modes (i.e. at an 0(1) distance 
downstream from the critical £ -location); in such cases the most dangerous mode has too large a 
growth rate and the wavenumber will not be close enough to a “neutral’ value for weakly-nonlinear 
theory to be immediately applicable. However, it can be argued that viscous spreading effects (or 
some other external effect) will reduce the growth rates to a size where a weakly nonlinear theory 
(based on unsteady critical-layer theory rather than viscous critical-layer theory) is appropriate. 
The papers by, for example, Michalke (1964); Crighton & Gaster (1976) and Hultgren (1992) 
support such an argument, which has been used in many recent papers concerned with flow 
stability eg. Goldstein & Leib (1988); Goldstein & Hultgren (1988); Goldstein & Leib (1989); 
Hultgren (1992); Wu, Lee & Cowley (1993). The evolution of the Hall-Horseman-modes for the 
non-marginal stability case is the subject of current study by the authors and will be reported on 
in due course. 

Our analysis shows that the disturbance amplitude satisfies the integro-differential equation 
(5.11/6.1). Experimental observations show that the linear growth of three-dimensional distur- 
bances to Gortler vortices is rapidly followed by the onset of turbulence. Such a scenario would 
be consistent with (6.1) if the positive sign were taken in that equation. The excellent agree- 
ment between the experimental measurements of Swearingen and Blackwelder (1987) and Hall 
and Horseman (1991) for the linear regime lead us to believe that this is indeed the case. However 
it is possible that the sign to be taken in (6.1,2) depends on the wavenumber and frequency of the 
marginally unstable mode and that some disturbances are destroyed by viscous effects. Therefore 
it is conceivable that the linearly growing disturbances are inhibited by viscosity for weak vortex 
states and grow explosively further downstream when the vortex state has been reinforced. 

Following the explosive growth of the disturbance new effects must come into play and viscosity 
will play a secondary role; see Wu, Lee and Cowley (1993) who were concerned with the growth 
of inviscid disturbances to Stokes layers. The extension of our work along the line followed by the 
latter and indeed other authors is made nontrivial by the fact that in our calculation the critical 
layer is not flat. 
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Figure Captions 


Figure 1. The growth rates of the two most dangerous odd modes at a position 100cm from the 
leading edge in the experiment of Swearingen and Blackwelder (1987). 

Figure 2. The contours of the streamwise and normal velocity components for the most dangerous 
odd mode 100 cm from the leading edge in the experiment of Swearingen and Blackwelder (1987). 

Figure 3. The various regions of the three-dimensional boundary-layer flow: Ia,b The outer 
flow; IIa,b The diffusion layers of relative thickness (^(ife -1 / 4 // -1 / 2 ); III Critical layer of relative 
thickness 0(Re ~ 1/,<3 ); IV Outer wall lay er of relative thickness 0(Re 1/,<3 ); V Inner wall layer of 
relative thickness 0(i?e -1 / 4 ). Note that (i), these thicknesses are given relative to the boundary- 
layer thickness 0(i?e -1 / 2 ), and (ii), the critical and diffusion layers are distinct from the wall being 
centered around the level Y = f(x, Z ). 
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